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COLD SCALAR-TENSOR BLACK HOLES: 
CAUSAL STRUCTURE, GEODESICS, STABILITY 

K.A. Bronnikov 1 , G. Clement 2 , CP. Constantinidis 3 and J.C. Fabris[] 

Departamento de Fisica, Universidade Federal do Espirito Santo, Vitoria, Espirito Santo, Brazil 

We study the structure and stability of the recently discussed spherically symmetric Brans-Dicke black-hole type 
solutions with an infinite horizon area and zero Hawking temperature, existing for negative values of the coupling 
constant uj . These solutions split into two classes: Bl, whose horizon is reached by an infalling particle in a finite 
proper time, and B2, for which this proper time is infinite. Class Bl metrics are shown to be extendable beyond 
the horizon only for discrete values of mass and scalar charge, depending on two integers m and n. In the case 
of even m — n the space-time is globally regular; for odd m the metric changes its signature as the horizon is 
crossed. Nevertheless, the Lorentzian nature of the space-time is preserved, and geodesies are smoothly continued 
across the horizon, but, when crossing the horizon, for odd m timelike geodesies become spacelike and vice versa. 
Problems with causality, arising in some cases, are discussed. Tidal forces are shown to grow infinitely near type Bl 
horizons. All vacuum static, spherically symmetric solutions of the Brans-Dicke theory with uj < —3/2 are found 
to be linearly stable against spherically symmetric perturbations. This result extends to the generic case of the 
Bergmann- Wagoner class of scalar-tensor theories of gravity with the coupling function u>{4>) < —3/2. 



1. Introduction 

In the recent years there has been a renewed inter- 
est in scalar-tensor theories (STT) of gravity as viable 
alternatives to general relativity (GR), mostly in con- 
nection with their possible role in the early Universe: 
they provide power-law instead of exponential infla- 
tion, leading to more plausible perturbation spectra 
Another aspect of interest in STT is the possi- 
ble existence of black holes (BHs) different from those 
well-known in GR. Thus, Campanelli and Lousto Q 
pointed out among the static, spherically symmetric 
solutions of the Brans-Dicke (BD) theory a subfam- 
ily possessing all BH properties, but (i) existing only 
for negative values of the coupling constant u> and (ii) 
with horizons of infinite area. These authors argued 
that large negative uj are compatible with modern ob- 
servations and that such BHs may be of astrophysical 
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relevance^. Indeed, the post-Newtonian parameters, 
determining the observational behaviour of STT in the 
weak field limit, depend crucially on the absolute value 
of uj rather than its sign [Q, || . 

In Ref. H it was shown, in the framework of a gen- 
eral (Bergmann- Wagoner) class of STT, that nontriv- 
ial BH solutions can exist for the coupling function 
oj(4>) + 3/2 < 0, and that only in exceptional cases 
these BHs have a finite horizon area. In the BD theory 
(uj = const) such BHs were indicated explicitly; they 
have infinite horizon areas and zero Hawking temper- 
ature ("cold BHs"), thus confirming the conclusions of 
These BHs in turn split into two subclasses: Bl, 
where horizons are attained by infalling particles in a 
finite proper time r, and B2, for which r is infinite. 
These results are briefly presented in Sections 2 and 3. 

The static region of a type B2 BH is geodesically 
complete since its horizon is infinitely remote and ac- 
tually forms a second spatial asymptotic. For type Bl 
BHs the global picture is more complex and is dis- 
cussed here in some detail (Sec. 4). It turns out that 
the horizon is generically singular due to violation of 
analyticity, despite the vanishing curvature invariants. 
Only a discrete set of Bl-solutions, parametrized by 

5 For brevity, we call BHs with infinite horizon areas type B 
BHs H, to distingish them from the conventional ones, with 
finite horizon areas, to be called type A. 



2 



K.A. Bronnikov, G. Clement, CP. Constantinidis and J.C Fabris 



two integers m and n , admits a Kruskal-like extension, 
and, depending on their parity, four different global 
structures are distinguished. Two of them, where m— n 
is even, are globally regular, in two others the region 
beyond the horizon contains a spacelike singularity. 

Subsec. 4.4 describes the behaviour of geodesies for 
different cases of BD BH metrics. It turns out that for 
odd m, when crossing the horizon, timelike geodesies 
become spacelike and vice versa, leading to problems 
with causality. Thus, in a family of BHs there ap- 
pear closed timelike curves, whose existence may be 
avoided by assuming a "helicoidal" analytic extension 
of the space-time manifold. Moreover, it is shown by 
considering the geodesic deviation equations that near 
a horizon with an infinite area all extended bodies are 
destroyed (stretched apart) by unbounded tidal forces. 

Sec. 5 discusses the stability of STT solutions un- 
der spherically symmetric perturbations. In this case 
the only dynamical degree of freedom is the scalar 
field, and the perturbation analysis reduces to a single 
wave equation whose radial part determines the system 
behaviour. Under reasonable boundary conditions, it 
turns out that the BD solutions with to < —3/2 are 
linearly stable, and this result extends to similar solu- 
tions of the general STT provided the scalar field does 
not create new singularities in the static domain. For 
the case u> > —3/2 the stability conclusion depends on 
the boundary condition at a naked singularity, which 
is hard to formulate unambiguously. 

We can conclude that, in STT with the anoma- 
lous sign of the scalar field energy, vacuum spherically 
symmetric solutions generically describe stable BH-likc 
configurations. Some of them, satisfying a "quantiza- 
tion" condition, are globally regular and have peculiar 
global structures. We also conclude that, due to infi- 
nite tidal forces, a horizon with an infinite area con- 
verts any infalling body into true elementary (point- 
like) particles, which afterwards become tachyons. 

A short preliminary account of this work has been 
given in [Q. In our next paper we intend to discuss 
similar solutions with an electric charge. 

2. Field equations 

The Lagrangian of the general (Bergmann- Wagoner) 
class of STT of gravity in four dimensions is 



4>R 



(i) 



where uj((f>) is an arbitrary function of the scalar field 
<p and L m is the Lagrangian of non-gravitational mat- 
ter. This formulation (the so-called Jordan confor- 
mal frame) is commonly considered to be fundamental 
since just in this frame the matter energy-momentum 



tensor T{f obeys the conventional conservation law 
V Q T^ = 0, giving the usual equations of motion (the 
so-called atomic system of measurements) . In particu- 
lar, free particles move along geodesies of the Jordan- 
frame metric. Therefore, in what folows we discuss the 
geometry, causal structure and geodesies in the Jordan 
frame. 

We consider only scalar- vacuum configurations and 
put L m = 0. 

The field equations are easier to deal with in the 
Einstein conformal frame, where the transformed scalar 
field ip is minimally coupled to gravity. Namely, the 
conformal mapping <7 M „ = <^ _1 g^ transforms Eq. ([!]) 
(up to a total divergence) to 



L 



-g[R 



e = sign(w + 3/2), — = 

The field equations are 



3/2 



1/2 







(2) 
(3) 

(4) 



where we have suppressed the bars marking the Ein- 
stein frame. The value e = +1 corresponds to nor- 
mal STT, with positive energy density in the Einstein 
frame; the choice s = — 1 is anomalous. The BD the- 
ory corresponds to the special case u> — const , so that 



<f> = exp^/Vk + 3 / 2 D- 



(5) 



Let us consider a spherically symmetric field sys- 
tem, with the metric 



rfs| = c 2 ~<dt 2 - . 

dn 2 = do 2 



>du 2 



^dn 2 



sin 2 6W$ 2 , 



(6) 



where E stands for the Einstein frame, u is the radial 
coordinate, a, [3, 7 and the field ip are functions of 
u and t . Preserving only linear terms with respect to 
time derivatives, we can write Eqs. (§) in the following 
form: 



e 2a R° = e 2a - 2 "*{a + 2'l3) 

_[ 7 " + 7 ' (7 '_ a ' +2 /3')] 
e 2a R\ = e 2a ~ 2 ~<& 



0: 



(7) 



- [7" + 2/3" + 7 ' 2 + 2(3' 2 - a\i + 213')] = W 2 ; (8) 



2a D 2 „2a-2/3 



e 



+ c 



2a-27 



ft 



-\P" + f3'{ 1 '-a' + 2[3')} = Q; (9) 
i?oi = 2[[3' + (3[3' - a(3' ~ Pi] =-e<ptf\ (10) 



f^gUtp = c-"< +a+2l3 <p - 1 c 



where dots and primes denote, respectively, d/dt and 
d/du. 

Up to the end of Sec. 4 we will restrict ourselves to 
static configurations. 
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3. Black holes in scalar-tensor theories 

The general static, spherically symmetric scalar- vacu- 
um solution of the theory (pi) is given by M a 



ds 2 



1 



-2bu 



Jlbu 



dt- 



s 2 (k, u) 



du 2 



s 2 (k, u) 



dW 



tp = Cu + (pa, 



C, ipo — const, 



(12) 
(13) 



where J denotes the Jordan frame, u is the harmonic 
radial coordinate in the static space-time in the Ein- 
stein frame, such that a(u) = 2(3(u) + r y(u) , and the 
function s(k,u) is 

sinh/cu, k > 0, 
u, k = 0, (14) 

fc _1 sin ku, k < 0, 

The constants b, k and C (the scalar charge) are re- 
lated by 



2k 2 sign k = 2b 2 + eC 2 



(15) 



The range of u is < u < w max , where u — cor- 
responds to spatial infinity, while w max may be finite 
or infinite depending on k and the behaviour of <p((p) 
described by (§). In normal STT (e = +1), by @, 
we have only k > , while in anomalous STT k can 
have either sign. 

According to the previous studies [Q, [|, all these 
solutions in normal STT have naked singularities, up 
to rare exceptions when the sphere u — oo is regular 
and admits an extension of the static coordinate chart. 
An example is a conformal scalar field in GR viewed 
as a special case of STT, leading to BHs with scalar 
charge || |). Even when it is the case, such configu- 
rations are unstable due to blowing-up of the effective 
gravitational coupling jl0| . 

In anomalous STT (e = — 1 ) the solution behaviour 
is more diverse and the following cases without naked 
singularities can be found: 

1 . k > . Possible event horizons have an infinite area 
(type B BHs), i.e. 522 ^00 as r — > 2k. In BD theory, 
after the coordinate transformation 



-2kt 



(16) 



e — = 1 - 2k/r = P(r) 
the solution takes the form 
ds] = P~^ds% 

= p-t (P a dt 2 - p- a dr 2 - P 1 - a r 2 dn 2 ' 
4> = I* (17) 
with the constants related by 

(2lu + 3)£ 2 = 1 - a 2 , a = b/k. (18) 



The allowed range of a and £, providing a horizon 
without a curvature singularity at r = 2k , is 



a > 1, 



a > £ > 2 



(19) 



(Eqs. are valid for uj > -3/2 as well, but then 

a < 1, leading to a naked singularity.) 

For £ < 1 particles can arrive at the horizon in a 
finite proper time and may eventually (if geodesies can 
be extended) cross it, entering the BH interior (type Bl 
BHs ||). When £ > 1, the sphere r = 2k is infinitely 
far and it takes an infinite proper time for a particle 
to reach it. Since in the same limit 322 - * 00, this 
configuration (a type B2 BH ||) resembles a wormhole. 

2. k — 0. Just as for k > 0, in a general STT, only 
type B black holes are possible ||), with 322 -> 00 as 
u — ► 00 . In particular, the BD solution is 

- 2bu 'du 2 
72 



ds 1 



-2bu 



dt 2 



■dri- 



ll* 



s 2 (w + 3/2) = -26 2 



(20) 



The allowed range of the integration constants is b > 
0, 2b > s > —2b. This range is again divided into 
two halves: for s > we deal with a type Bl BH, for 
s < with that of type B2 ( s = is excluded since it 
leads to GR). 

3. k < . In the general STT the metric ( |l2| ) describes 
a wormhole, with two flat asymptotics at u = and 
u = 7r/|fc|, provided <f> is regular between them. In 
exceptional cases the sphere u max = 7r/|fc| may be an 
event horizon, namely, if <j> ~ 1/Au 2 , Au = \u— M max | • 
In this case it has a finite area (a type A black hole) 
and 



uj((f>) +3/2 



-0 



u — > u„ 



(21) 



The behaviour of goo and gn near the horizon is then 
similar to that in the extreme Reissner-Nordstrom so- 
lution. 

In the BD theory we have only the wormhole solu- 
tion 



ds = c 



-21, 



u dt 2 - 



k 2 e 2bu ( k 2 du 2 



sin 2 ku V sin 2 ku 



dfl' 



s 2 {uj + 3/2) = -k 2 -2b 2 



(22) 



with masses of different signs at the two asymptotics. 

For all the BH solutions mentioned, the Hawking 
temperature is zero. Their infinite horizon areas may 
mean that their entropy is also infinite; however, a 
straightforward application of the proportionality re- 
lation from GR between entropy and horizon area, is 
here hardly justified; a calculation of BH entropy is a 
separate problem, discussed in a large number of recent 
works from various standpoints. 
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4. Analytic extension and causal 
structure of Brans-Dicke black 
holes 

4.1. Extension 

Let us discuss possible Kruskal-like extensions of type 
Bl BH metrics @ (fc > 0) and @ (k = Q) of the 
BD theory. 

For (|i~7|), with a>l>£>2 — a, we introduce, as 
usual, the null coordinates v and w : 



t + x. 



w 



t 



dcf 



p- a dr (23) 



where x — > oo as r — > oo and a; — > — oo as r — > 2fc. 
The asymptotic behaviour of x as r — > 2fc (P — > 0) 
is i a — P 1_a , and in a finite neighbourhood of the 
horizon P = one can write 



a;=i(t;-«;) = -iP 1 -/(P), 



(24) 



where f(P) is an analytic function of P, with /(0) = 
4fc/(a- 1): 



/(P) = -4fc]T^±^r<. 



<?=0 



(? - a + 1 



(25) 



Then, let us define new null coordinates V < and 
VF > related to v and w by 



n = const. (26) 



The mixed coordinate patch (V, w) is defined for v < 
(t < — jc) and covers the whole past horizon v = — oo. 
Similarly, the patch (v, W) is defined for w > (i > 
x) and covers the whole future horizon w = +oo. So 
these patches can be used to extend the metric through 
one or the other horizon. 

Consider the future horizon. As is easily verified, 
a finite value of the metric coefficient g v w at W = 
is achieved if we take n + 1 = (a — 1)/(1 — £) , which 
is positive for a > 1 > £. In a finite neighbourhood 
of the horizon, the metric ( |l7|) can be written in the 
coordinates (u, 14 7 ) as follows: 

= -(n+ 1)/^tt . (i _ u W rn+1 )-^Tdw(iW 



4fc 2 



f-^tx . (l - w Ty n+1 )^TT^" m dfi 



(1-P)2 



where 



m= (o-l + 0/(l -0- 



(27) 



(28) 



The metric (p7| ) can be extended at fixed v from W > 
to T4 7 < only if the numbers n + 1 and m are both 




a; — > — oc 




X — > — oo 



Figure 1: Extensions through the future horizon for odd 
(a) and even (b) values of n . Thick vertical and horizontal 
arrows show the growth direction of the x coordinate in 
the corresponding regions. 



integers (since otherwise the fractional powers of neg- 
ative numbers violate the analyticity in the immediate 
neighbourhood of the horizon) . This leads to a discrete 
set of values of the integration constants a and £ : 

m + 1 to — n — 1 , . 

a= , £= . 29 

to — n m — n 

where, according to the regularity conditions (|i9|), 

to > n > . Excluding the Schwarzschild case to = 

n + 1 (£ = 0), which corresponds from (|22[) to a = 1 

(to = 0), we see that regular BD BHs correspond to 

integers to and n such that 



2 > n > 0. 



(30) 



An extension through the past horizon can be per- 
formed in the coordinates (V, w) in a similar way and 
with the same results. 

It follows that, although the curvature scalars van- 
ish on the Killing horizon P = 0, the metric cannot 
be extended beyond it unless the constants a and £ 
obey the "quantization condition" (^9|), and is generi- 
cally singular. The Killing horizon, which is at a finite 
affine distance, is part of the boundary of the space- 
time, i.e. geodesies and other possible trajectories ter- 
minate there. Similar properties were obtained in a 
(2+l)-dimensional model with exact power-law met- 
ric functions fllj and in the case of black p-branes 

We have thus obtained a discrete family of BH so- 
lutions whose parameters depend on the two integers 
to and n. The corresponding parameters describing 
the asymptotic form of the solution, the active gravi- 
tational mass M and the scalar charge S , defined by 

1 2GM fl\ , 1 S (Y 

.goo = 1 f- o| - K = 1 
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where G is the Newtonian gravitational constant, are 
"quantized" according to the relations 



GM 



-2k- 



m — n — 1 



(31) 



The constant k , specifying the length scale of the so- 
lution, remains arbitrary. On the other hand, the cou- 
pling constant co takes, according to (flcf), discrete val- 
ues: 



|2w + 3| 



(2m -n+ l)(n + 1) 
(m — n — l) 2 



(32) 



The k = solution (|2C|) of the BD theory also has 
a Killing horizon (it — > oo) at finite geodesic distance 
if s > 0. However, this space-time does not admit a 
Kruskal-like extension and so is singular. The reason 
is that in this case the relation giving the tortoise-like 
coordinate x, 



g26u e 2 ^ u 

— du =2b^ F{u) 



(33) 



(where F(u) is some function such that F(oo) = 1) 
cannot be inverted near u = oo to obtain u as an 
analytic function of x. 

4.2. Geometry and causal structure 

To study the geometry beyond the horizons of the met- 
ric ( |l7|) , or ( p7| ) , let us define the new radial coordinate 
P by 



P EE e- 2kU EE 1 — — — EE p n 



(34) 



The resulting solution (P~T|) , defined in the static region 
I (P > 0, p > 0), is 



ds 2 



n+2 M 2 _ 4k 2 (m-n) 2 

4 " 



p- 1 ~ df 

4fc 2 
(1-Py 



(1-P)< 



"- 2 rfp 2 

m—n— 1 



(35) 



By (|24|), p is related to the mixed null coordinates 
(v, W) by 

p(v, W) = W [/(P)] 1/(n+1) [l - vW n+1 ]- 1/[n+1) - 

(36) 

This relation and a similar one giving p(V, w) show 
that when the future (past) horizon is crossed, p varies 
smoothly, changing its sign with W (V). For p < 
the metric (|35| ) describes the space-time regions be- 
yond the horizons. 

To construct the corresponding Penrose diagrams, 
it is helpful to notice that by (|24|) the radial coordinate 
x is related to p by 



(37) 



so that for odd n the horizon as seen from region II 
(p < 0) also corresponds to x — > — oo. On the other 
hand, in the 2-dimensional metric ds\ = p n+2 (dt 2 — 
dx 2 ) , for p < the coordinate x is timelike, hence in 
region II beyond the future horizon (with respect to 
the original region I) x — > — oo means "down". A new 
horizon for region II joins the picture at point O - 
see Fig. 1(a). For even n, when x in region II remains 
a spatial coordinate and the coordinate ( |37| ) changes 
its sign when crossing the horizon, a new horizon joins 
the old one at the future infinity point of region I — 
Fig. 1(b). These considerations are easily verified by 
introducing null coordinates in region II, similar to 
v and w previously used in region I. Continuations 
through the past horizons are performed in a similar 
manner. The resulting causal structures depend on the 
parities of m and n . 

la. Both m and n are even, so P{p) is an even func- 
tion. The two regions p > and p < are isometric 
{9^v(—p) — 9nv{p))i an d the Penrose diagram is simi- 
lar to that for the extreme Kerr space-time, an infinite 
tower of alternating regions I and II (Fig. 2, left). All 
points of the diagram, except the boundary and the 
horizons, correspond to usual 2-spheres. 

lb. Both m and n are odd; again P(p) is an even 
function, but regions I and II are now anti-isometric 
i9nu(~p) = ~9iJ,v{p))- The metric tensor in region 

II (p < 0) has the signature ( h ++) instead of 

(H ) in region I. Nevertheless, the Lorentzian 

nature of the space-time is preserved. The Penrose 
diagram is shown in Fig. 3, left. In both cases la and 
lb the maximally extended space-times are globally 
regular^. 

2. rn — n is odd, i.e. P(p) is an odd function; more- 
over, P — > — oo (r — > 0) as p — > — oo, so that the 
metric ([55]) is singular on the line p = — oo, which 
is spacelike. The resulting Penrose diagrams are sim- 
ilar to those of the Schwarzschild space-time in case 
2a (n odd, m even), and of the extreme (e 2 = to 2 ) 
Reissner-Nordstrom space-time in case 2b ( n even, to 
odd, Fig. 4). In the latter subcase, however, the 4- 
dimensional metric changes its signature when crossing 
the horizon, similarly to case lb. 

4.3. Type B2 structure 

Let us briefly consider the case B2: fc>0, a > £ > 1. 
As before, the metric is transformed according to ( ^3|) 
( |26] ) and at the future null limit (now infinity rather 
than a horizon, therefore we avoid the term "black 



e A globally regular extension of an extreme dilatonic black 
hole, with the same Penrose diagram as in our case la, was 



discussed in |12| . 
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Figure 2: The Penrose diagram and the effective potential for geodesies for a BH with m and n both even. Curve 1 
is V v (p) , curve 2 is Vl(p) and curve 3 is V(p) for a nonradial timelike geodesic. The curves and "energy levels" E1-E5 
correspond to different kinds of geodesies as described in Subsec. 4.4, item la. 



hole" ) we again arrive at ( p7| ) , where now W — > oo as 
P — > . The asymptotic form of the metric as W — > oo 
is 

ds 2 = -ddv dW - C 2 W m dn 2 (38) 

where C± t 2 are some positive constants, while the con- 
stant to, defined in (|2g|), is now negative. A further 
application of the v -transformation (|2^) at the same 
asymptotic, valid for any finite v < , leads to 

ds 2 = -ai-vfr-w-VdVdW - c 2 w- m dn 2 . 

(39) 

If we now introduce new radial (R) and time (T) coor- 
dinates by T = V + W and R = T — W , in a spacelike 
section T = const the limit R — > — oo corresponds 
to simultaneously V — > — oo and W — * +oo , with 
|V| ~ W, and the metric ( |39| ) turns into 

da 2 = 4Ci(-i?) (Q -« )/( «- 1) (dT 2 -di? 2 ) 

- C 2 (-ii)- ro dfi 2 . (40) 

Evidently, this asymptotic is a nonflat spatial infin- 
ity, with infinitely growing coordinate spheres and also 
infinitely growing goo j i-e. , this infinity repels test par- 
ticles. 

A Penrose diagram of a B2 type configuration co- 
incides with a single region I in any of the diagrams; 
all its sides depict null infinities, its right corner corre- 
sponds to the usual spatial infinity and its left corner 
to the unusual one, described by the metric j4C|). The 
latter has been obtained here by "moving along" the 



future null infinity W — > oo , but the same is evidently 
achieved starting from the past side. 

4.4. Geodesies 

Let us now return to type Bl BHs and study test parti- 
cle motion in their backgrouds, described by geodesies 
equations. 

One can verify that all geodesies are continued 
smoothly across the horizons, even in cases lb and 2b 
when the metric changes its signature (the geodesic 
equation depends only on the Christoffel symbols and 
is invariant under the anti-isometry -— > —g^v ) . 

We will use the metric (|35|). Then the integrated 
geodesic equation for arbitrary motion in the plane 9 = 
tt/2 reads: 

Ak [t~ P f ? + W " +2 + S (1 - P)2pm+n+2 = ^ 

(41) 

where p = dp/dX, A being an affine parameter such 
that ds 2 = r/dX 2 , with 77 = +1, 0, —1 for timelike, 
null or spacelike geodesies; E 2 and L 2 are constants 
of motion associated with the timelike and azimuthal 
Killing vectors and, correspondingly, with the particle 
total energy and angular momentum. Eq. ( f4"l| ) has the 
form of an energy balance equation, with the effective 
potential 

V(p)=V v (p) + V L (p) (42) 
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Figure 3: The Penrose diagram and the effective potential for geodesies for a BH with m and n both odd. Curve 1 is 
V v (p) , curve 2 is Vl(p) and curve 3 is V(p) for a nonradial timelike geodesic. The curves and "energy levels" E1-E7 
correspond to different kinds of geodesies as described in Subsec. 4.4, item lb. 



where Vq and Vr, are the respective terms in the left- 
hand side. 

One should note that although the coordinate p 
belongs to the static frame of reference, one can use it 
(and consequently Eq. ([ll])) to describe geodesies that 
cross the horizon since in a close neighbourhood of the 
horizon p = it coincides (up to a positive constant 
factor) with a manifestly well-behaved coordinate V 
or W and, on the other hand, Eq. ( f4l| ) reads simply 
p 2 = const • E 2 and is thus also well-behaved. 

Let us discuss the four possible cases according to 
the parity of m and n and the corresponding signature 
of the metric ([55]) for p < : 

la: m even, n even, (H ). The range of p is 

(— 1,+1). The coefficient (1 — P)^ 1 of the kinetic term 
and the potential in (^lj) are both symmetric under 
the exchange p — » —p. The potential V(p) is shown 
in Fig. 2: curve 1 depicts V n (p) for 77 = 1, i.e., the 
potential for radial timelike geodesies; curve 2 shows 
Vl , the angular momentum dependent part of V(p) , 
and curve 3 shows their sum for certain generic val- 
ues of the motion parameters. Depending on the value 
of E, geodesic motion can be symmetrical with suc- 
cessive horizon crossings from one region to the next 
isometrical region without reaching past or future null 
infinity (El: see "energy levels" in Fig. 2, right and the 
corresponding curves in Fig. 2, left), or starting from 
a past timelike infinity in region I and reaching a fu- 
ture timelike infinity in region II (E4). Some nonradial 
timclike (E2, E3) and null (E5) geodesies remain in a 



single region, corresponding to bound (E2) or unbound 
(E3) particle orbits near the BH or photons passing it 
by (E5). The existence of nonradial trajectories like 
El for any value of L is here connected with the infi- 
nite value of e' 3 at the horizon, creating a minimum 
of Vl- Radial null geodesies (V(p) = 0) correspond, 
as always, to straight lines tilted by 45° (unshown). 

lb: m odd, n odd, ( — h + +). The one-dimensional 
dependence p(A) is qualitatively the same for null 
geodesies (77 = 0). However, the global picture is 
drastically different, see Fig. 3. In particular, type 
El null geodesies which periodically cross the horizon 
necessarily go repeatedly through all four regions of 
the Penrose diagram, so that their projections onto a 
2-dimensional plane 9 = const, $ = const are closed; 
they will be even closed in the full 4-dimensional space- 
time for some discrete values of the angular momen- 
tum L 2 . Thus such BH space-times contain closed 
null geodesies, leading to causality violation. 

However, this problem can be avoided if we choose, 
instead of the simplest, one-sheet maximal analytic ex- 
tension corresponding to the planar Penrose diagram 
of Fig. 3, left, a "helicoidal" analytic extension con- 
structed in the following way: starting from a given re- 
gion I and proceeding with the extension counterclock- 
wise, after 4 steps we come again to a region I isometric 
to the original one, but do not identify these mutually 
isometric regions and repeat the process indefinitely. 
The same process is performed in the clockwise direc- 
tion. The resulting Penrose diagram is a Riemann sur- 
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Figure 4: The Penrose diagram and the effective potential for geodesies for a BH with n even and m odd. Curve 1 is 
V v (p) , curve 2 is Vl(p) and curve 3 is V(p) for a nonradial timelike geodesic. The curves and "energy levels" E1-E7 
correspond to different kinds of geodesies as described in Subsec. 4.4, item 2b. 



face with a countable infinity of sheets similar to that in 
Fig. 3, left, cut along one of the horizons. Such a struc- 
ture no longer exhibits causality violation^. In general, 
when crossing a horizon, a null geodesic remains null, 
but timelike geodesies become spacelike and vice versa, 
since the coefficient r\ in Eq. , being an integral of 
motion for a given geodesic, changes its meaning in 
transitions between regions I and II: a geodesic with 
T] = 1 is timelike in region I and spacelike in region II, 
and conversely for rj = — 1 . 

For nonnull geodesies the potential is now asym- 
metric: thus, for trajectories which are timelike in 
region I (rj = +1), it becomes attractive for p < 
(Fig. 3, the notations coincide with those of Fig. 2). 
These geodesies become spacelike in region II and 
generically extend to spacelike infinity as A — > oo 
(E4 in Fig. 3). This is true for all radial geodesies and 
part of nonradial ones; however, nonradial geodesies 
with small E (near the minimum of curve 3 at p = ) 
are of the type El, quite similar to null El trajec- 
tories that have been just discussed. A new type of 
tachyonic motion as compared with item la is E6 
shown in Fig. 3. There are also circular unstable 
geodesies with i] — +1 in region II, with t = const 
and p = pa = — [m/(3m — 2n)] 1 ^" 1- "^ (points in the 



7 Strictly speaking, such a process might be applied to the 
Schwarzschild and Rindler space-times, with possible identifica- 
tion of isometric regions after a finite number of steps. This is, 
however, unnecessary since there a causality problem like ours 
does not exist. 



Penrose diagram), such that V(po) — V'(pa) = 0, 
E = 0. All these spacelike geodesies have full ana- 
logues with r\ = — 1 in region I. 

The whole space-time possesses full symmetry un- 
der an exchange between regions I and II, correspond- 
ing to rotations of the Penrose diagram of Fig. 3 by 
odd (in addition to even) multiples of the right angle, 
accompanied by a change of sign of r\ so that geodesies 
keep their timelike or spacelike nature. 

2a: m even, n odd, ( 1 ). The range of p is 

now (— oo,+l). Both parts of the effective potential 
become negative and monotonic at p < 0, so that all 
geodesies entering the horizon terminate at the space- 
like singularity p — — oo, as in the Schwarzschild case. 
Thus the whole qualitative picture of test particle mo- 
tion, as well as the Penrose diagram, coincide with the 
Schwarzschild ones. 

2b: m odd, n even, (H h +) . The potential V v 

[rj — 1 ) is positive-definite (as in the extreme Reissner- 
Nordstrom case), as shown in Fig. 4, so that all radial 
timclike geodesies avoid the singularity, crossing the 
horizon either twice (E6), or indefinitely (El). For non- 
radial motion, the summed potential V(p) — > — oo as 
p — > — oo , whatever small is I? , therefore all geodesies 
that are timelike in region I (and become spacelike in 
region II) reach the spacelike singularity^ p — — oo 

8 Due to the metric signature change, in the Penrose diagram 
of Fig. 4, just as in Fig. 3 for the case lb, the time direction is 
vertical in regions I and horizontal in regions II. 
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(levels and curves E4 and E7 in Fig. 4) , except those 
with small E depicted as El. 

One can conclude that the unusual nature of metric 
of Bl type BHs creates some unusual types of particle 
motion. Some of them even exhibit evident causal- 
ity violation — such as an observer receiving messages 
from his or her own future — which can be avoided by 
assuming a more complicated space-time structure. 

Another paradox, also related to causality, arises if 
we follow in cases lb or 2b the fate of a hypothetical 
(timelike) observer who has crossed the horizon and 
finds him(her)self in a region II where his (her) proper 
time is now spacelike as viewed by a resident observer 
(whose timelike geodesic is entirely contained in region 
II), and can be reversed by a simple coordinate change. 

However, one may suspect that, the horizon area 
being infinite, any extended body, and an observer 
in particular, will have been infinitely stretched apart 
and destroyed before actually crossing the horizon. To 
check this, consider for instance a freely falling ob- 
server whose center of mass follows a radial geodesic 
in the plane = tt/2. The separation n a between 
this geodesic and a neighbouring radial geodesic varies 
according to the law of geodesic deviation 



D n° 
U 2 ~ 



R a 0lS u^n i u & = 0. 



(43) 



For the four-velocity of the center of mass we have 
u° = Eg 00 and u°uq + u x u\ = 1, so that near the 
horizon u» ~ {Eg 00 , E{-g m g 11 ) 1 ' 2 , 0, 0) . We obtain 
for the relative azimuthal acceleration near the horizon 



1 D 



2„3 



n 3 d\ 2 



= {R 30 30 u°u + tfWui) 



-E 2 g 00 (R 30 - 



30 



i? 31 31 ) = E Z R"/R, (44) 



and a similar equation for the deviation n 2 in the 9 
direction. Here R 2 = | (?22 1 , R" — d 2 R/dp 2 , and p is 
a radial coordinate such that the Jordan-frame metric 
functions are related by <?oo5n = const; this condition 
is valid near the horizon for our coordinate p defined 
in (g). 

The azimuthal geodesic deviation (|44|) which, due 
to its structure, is insensitive to radial boosts and is 
thus equally applicable to the static frame of reference 
and to the one comoving with the infalling body, agrees 
with similar relations given by Horowitz and Ross JH| . 
In the case of the Schwarzschild metric, R" /R — and 
the tidal force (given by the terms that we have ne- 
glected) is finite. In the case of the examples discussed 
in R"/R is negative and large, i.e. geodesies con- 
verge and physical bodies are crushed as they approach 
the horizon, as by a naked singularity, hence the name 
"naked black holes" given to these spacetimes in 
On the contrary, in the case of the cold black hole 



metric (p5|), R" / R — > +oo (as p~ 2 ), i.e. geodesies di- 
verge; the resulting infinite tidal forces pull apart all 
extended objects, e.g. any kind of clock, approaching 
the horizon. Only true elementary (pointlike) parti- 
cles, resulting from the destruction of the falling body, 
cross such a horizon to become tachyons[]. 



5. Stability 

Let us now study small (linear) spherically symmetric 
perturbations of the above static solutions (or static 
regions of the BHs), i.e. consider, instead of (p(u) , 



ip{u, t) — (p(u) + 6(p(u, t) 



(45) 



and similarly for the metric functions a, (3, 7, where 
<p(u) , etc., are taken from the static solutions of Sec. 2. 
We are working in the Einstein conformal frame and 
use Eqs. (|J) . The consideration applies to the whole 
class of STT ([[]); its different members can differ in 
boundary conditions, to be discussed below. 

In perturbation analysis there is the so-called gauge 
freedom, i.e. that of choosing the frame of reference 
and the coordinates of the perturbed space-time. The 
most frequently used frame for studying radial pertur- 
bations has been that characterized by the technically 
convenient condition 5/3 = Jul |l4|, It was ap- 
plied, however, to background configurations where the 
area function e" was monotonic in the whole range of 
u , or was itself used as a coordinate in the static space- 
time. Unlike that, in our study the configurations of 
utmost interest are type B black holes and wormholes, 
i.e. those having a minimum of e' 3 (a throat) at some 
value of u and infinite at both ends of the u range. 
At the throats, the equality 5(3 = is no longer a co- 
ordinate condition, but a physical restriction, forcing 
the throat to be at rest. It can be explicitly shown 
that the condition 5(3 = creates poles in the effective 
potential for perturbations, leading to their separate 
existence at different sides of the throat, i.e. the latter 
behaves like a wall. 

For these reasons, we have to use another gauge, 
and we choose it in the form 



5a = 25(3 + 5j, 



(46) 



extending to perturbations the harmonic coordinate 
condition of the static system. In this and only in this 
case the scalar equation ( |li"| ) for 5(p decouples from the 
other perturbation equations and reads 



] 5^-5ip" = 0. 



(47) 



9 Our cold black holes are thus counterexamples to the claim 
made in Jl3| ] that a smooth extension is not possible when tidal 
forces diverge on the horizon. 
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Since the scalar field is the only dynamical degree of 
freedom, this equation can be used as the master one, 
while others only express the metric variables in terms 
of dip , provided the whole set of field equations is con- 
sistent. That it is indeed the case, can be verified 
directly. Indeed, under the condition (fl6|), the four 
equations (@)-([l0|) for perturbations take the form 



(AS/3 + 5jJ- 5i' = 0; 



(48) 



e Af3 {25[3 + 8 1 y-25{3" -8^" 

- 4(0' - j')S(3' + 05*/' = 2eip'8ip'; (49) 

- 5(3" + 2 e 2/3+27 (6/3 + &y) = 0; (50) 
p' -f3' [5$ + 5^-^50 = -\eip' Sip, (51) 

where a, (3, 7, tp satisfy the static field equations. 

Eq. (^lj) may be integrated in t and further differ- 
entiated in u; the result turns out to be proportional 
to a linear combination of the remaining Einstein equa- 
tions ([l8|)-(|5(]). On the other hand, the quantities Sip 
and Sip" can be calculated from (^8|)-([5l]), resulting in 
p?|). Therefore we have three independent equations 
for the three functions Sip , S f3 and £7 . 

The following stability analysis rests on Eq. ([47j). 
The static nature of the background solution makes it 
possible to separate the variables, 



Sip = ip(u) e 1 ' 



(52) 



and to reduce the stability problem to a boundary- 
value problem for ip(u) . Namely, if there exists a non- 
trivial solution to d47| ) with lu 2 < 0, satisfying some 
physically reasonable boundary conditions at the ends 
of the range of u, then the static background system 
is unstable since perturbations can exponentially grow 
with t . Otherwise it is stable in the linear approxima- 
tion. 

Suppose — u! 2 = il 2 , £1 > . In what follows we use 
two forms of the radial equation ( ^7|) : the one directly 
following from (l52|), 



V/' - n 2 e 4/3 ("V = 0, 



(53) 



and the normal Liouville (Schrodinger-like) form 

d 2 y/dx 2 - [fl 2 + V(x)]y(x) = 0, 

V{x) = e- 4 ?{f3" - f3' 2 ). (54) 

obtained from (p3[) by the transformation 



ip(u) = y(x) e p , x 



(55) 



Here, as before, a prime denotes d/du. 

The boundary condition at spatial infinity ( u — > , 
x ~ 1/u — > +00) is evident: Sip — > 0, or ip — » 0. For 
our metric (|lj) the effective potential V(x) has the 
asymptotic form 



V(x) m 2b/x 3 , 



as 



+00, 



(56) 



hence the general solutions to (p4) and (B3J) have the 



asymptotic form 

Mx 1 



y 



Ci e 



c 2 e 



-fix 



{x — > +00), 



(57) 



ip - u( Cl e /" + c 2 e- 0/u ) (u^0), (58) 

with arbitrary constants c\ , c 2 . Our boundary condi- 
tion leads to cy = 0. 

For u — > ii max , where in many cases the back- 
ground field ip tends to infinity, a formulation of the 
boundary condition is not so evident. Refs. |lj, |l5| 
and others, dealing with minimally coupled or dilatonic 
scalar fields, used the minimal requirement providing 
the validity of the perturbation scheme in the Einstein 
frame: 



\5ip/ip\ < 00. 



(59) 



In STT, where Jordan-frame and Einstein-frame met- 
rics are related by g 3 = (l/0)gJL, it seems reasonable 
to require that the perturbed conformal factor l/4> be- 
have no worse than the unperturbed one, i.e. 



\84>/4>\ < 00. 



(60) 



An explicit form of this requirement depends on the 
specific STT and can differ from (|5^), for example, 
in the BD theory, where </> and ip are connected by 
(||), the requirement ( |60| ) leads to \8ip\ < 00. We will 
refer to (59) and ( |60| ) as to the "weak" and "strong" 
boundary condition, respectively. For configurations 
with regular </> and ip at u — > u max these conditions 
both give \Sip\ < 00. 

Let us now discuss different cases of the STT solu- 
tions. We will suppose that the scalar field <fi is regular 
for < u < it max , so that the conformal factor <$r 1 in 
( |l2] ) does not affect the range of the u coordinate. 

1. e = +1, k > 0. This is the singular solution of 
normal STT. As u — * +00 , f3 ~ (b — k)u — > — 00 , 
so that x tends to a finite limit and with no loss of 
generality one can put x — > . The effective potential 
V(x) then has a negative double pole, V ~ — 1/(4.t 2 ), 
and the asymptotic form of the general solution to ( |54"| ) 
leads to 



tp(u) w y(x)/\fx « (03 + c 4 lna;) (x — > 0). 



(61) 



The weak boundary condition leads to the require- 
ment \8ip/ip\ ~ |y|/ (-y/a^l l n ^|) < 00 , met by the general 
solution (|6l]) and consequently by its special solution 
that joins the allowed case (ci = 0) of the solution 
( |57j ) at the spatial asymptotic. We then conclude that 
the static field configuration is unstable, in agreement 
with the previous work Q . 

As for the strong boundary condition (|60|), proba- 
bly more appropriate in STT, its explicit form varies 
from theory to theory, therefore a general conclusion 
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cannot be made. In the special case of the BD theory 
the condition ( |60| ) means < oo as u — > +00. Such 
an asymptotic behavior is forbidden by (Eq. (|5^), ac- 
cording to which ip" '/ip > 0, i.e. the perturbation i/)(u) 
is convex and so cannot be bounded as u — > 00 for an 
initial value i/>(0) = (ci =0). We thus conclude that 
the static system is stable. 

We see that in this singular case the particular 
choice of a boundary condition is crucial for the stabil- 
ity conclusion. In the case of GR with a minimally cou- 
pled scalar field (lj] there is no reason to "strengthen" 
the weak condition that leads to the instability. In the 
BD case the strong condition seems more reasonable, 
so that the BD singular solution seems stable. For any 
other STT the situation must be considered separately. 

2. e = — 1, £; > 0. This case includes some singular 
solutions and cold black holes, as is exemplified for the 
BD theory in ©-(H). 

As u — > +00 , (3 ~ (b — k)u — > +00 , so that x — > 
—00 and V(x) w — l/(4a; 2 ) — > 0. The general solu- 
tion to Eq. (|54|) again has the asymptotic form ( |57| ) for 
x — ► —00 . The weak condition (|59| ) leads, as in the pre- 
vious case, to the requirement |j/|/(y / jx[ln |x|) < 00, 
and, applied to (p7|), to C2 = 0. This means that the 
function ip must tend to zero for both u — > and 
u — ► 00 , which is impossible due to tp"/ip > 0. Thus 
the static system is stable. Obviously the more restric- 
tive strong condition ( |60| ) can only lead to the same 
conclusion. 

3. £ = —1, k = 0. There are again singular solutions 
and cold black holes. As u — > +00, x — > — 00 and 
again the potential V(x) — 1/4x 2 — * 0, leading to 
the same conclusion as in case 2. 

4. e = — 1, £; < 0. In the generic case the solution 
describes a wormhole, and in the exceptional case ( pl| ) 
there is a cold black hole with a finite horizon area. In 
all such cases, as u — > w ma x = ^/l^l 1 one nas x ~ * — 00 
and V ~ 1/|^| 3 — > 0, so that the stability is concluded 
just as in cases 2 and 3. 

Thus, generically, scalar- vacuum spherically sym- 
metric solution of anomalous STT are linearly sta- 
ble against spherically symmetric perturbations. Ex- 
cluded are only the cases when the field <fi behaves in 
a singular way inside the coordinate range < u < 
u max ; such cases should be studied individually. 
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